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INTRODUCTION TO THEMATIC MAPPER INVESTIGATIONS 


SECTION I: RADIOMETRY 

This section of the Introduction to Thematic Mapper (TM) 
Investigations Is an overview of papers which deal with radiometric 
characterization of the TM sensor. Spectral characteristics are summarized In 
the paper by Markham and Barker. Papers dealing primarily with geometric 
characteristics of TM are overviewed In Section 2 of this Introduction by 
Brian Markham. 

References In this Introduction will be In one of four forms: 

t P: Author (Example: P: Bernstein) 

"P" refers to a paper In one of the two TM volumes In these 
February, 1983 Proceedings, and the name of the author Is 
the first author. If there Is more than one paper by the 
same first author, the name of the second author will also 
be Included. 

• S:Author, p (Example: S:Barker, Ball et a1.«Vo1. 1.130-139) 

*5“ refers to a summary in one of the two volumes or 
abstracts In NASA Conference Publication 2326 (1984) 
entitled "Landsat-4 Science Investigations Sunmary 
Including December 1983 Workshop Results." The page 
location In the publication Is given by "p". In almost all 
cases there will be a reference to either the summary or to 
the papers In the proceedings, but not to both. 

• I: Author (Example: I:Mal11a) 

"I" refers to one of the papers In the May, 1984 Issue of 
IEEE Transactions on Geoscience and Remote Sensing entitled 
"Special Issue on Landsat-4." The table of contents of 
this Issue Is given at the end of this report. 


• Author, Year (Example: Engel, 1980) 

References to other papers are given by first author and 
year of publication, with the complete citation being at 
the end of this report. There Is no Intention to provide 
more than a few such references to other papers since this 
Is a report on the radiometric and geometric activities of 
the Landsat Image Data Quality Assessment (LIDQA) 
Investigators (Table 1). 

For most users, the primary scientific requirement Is to convert 
digital numbers (DN) on computer compatible tapes (CCT) to radiance or 
spectral radiance units. Minimum and maximum spectral radiances for the six 
reflective bands are given for radlometrlcally calibrated tapes for two 
periods of time: 

Scrounge-Era prior to August 1983 (S:Barker and Ball, Vol. I, 

130-139, Table 3) 

TIPS-Era after 15 January 1984 (SrBarker, Vol. I,^« 140-180, Table 

1; and P:Barker, Appendix 9.1 ). 

For both the reflective bands and the emissive thermal band, TM6, the 
actual dynamic range used In calibrating the digital Imagery Is contained In 
one of the fields of the CCT. For perfectly calibrated noise-free data, these 
post-calibration dynamic ranges are all the radiometric Information needed. 

Investigations reported here, and still In progress, are 
characterizing and monitoring the radiometric performance of the TM sensor on 
both Landsat-4 (TM/PF) and on Landsat-5 (TM/F). These efforts Include 
pre-launch calibration and In-orbit calibration of sensors. There has not yet 
been an opportunity to characterize the performance or optimize the parameters 
used In the radiometric ground processing of TM Imagery. This section Is 
divided Into three topics: characterizing pre-launch sensor performance, 
characterizing In-orbit sensor performance, and characterizing calibrated 
digital Image data quality. 

Pre-Launch Sensor Performance 


Summary charts on the TM sensor description and characteristics are 
given here (S:Enge1, Vol. 1, 41-61) and In a previous paper (Engel and 
Weinstein, 1983). Details of the absolute radiometric calibration using an 
Integrating sphere are given for l.andsat-4 TM/PF (P: Barker and Ball) along 
with summary results on gains and offsets for Landsat-5 TM/PF (P: Barker, 
Appendices 9.3 and 9.4). Many special tests were performed before launch, 
especially on Landsat-5 TM/F, and tapes from these tests provide an option for 
future re-cal 1b ration based on a knowledge of how to reduce certain systematic 
sources of variability that were not characterized until after launch. 

Re-cal Ibratlon of the Integrating sphere may also result In a further updating 
of the post-calibration dynamic range In the parameter file of the TM Image 
Processing System (TIPS). Independent measurement of absolute radiometric 
calibration Is being attempted by comparing In-orbit Imagery with simultaneous 
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observations from the ground (S:Cast1e, Vol. II, 15-19; S:Hov1s, Vol. I, 
181-185). Absolute calibration had a specified accuracy requirement of 101 of 
full scale In each band. 

In-Orbit Sensor Performance 


As part of the Investigations program, several scientific teams have 
had access to raw uncallbrated TM digital Imagery on engineering tapes. 
Including foreign ground stations with direct access to TM transmission from 
Landsat. Several radiometric studies have used uncallbrated digital Imagery, 
either CCT-BT tapes produced during the Scrounge-era or "unity" CCT-AT tapes 
from TIPS, In conjunction with the background shutter and Internal calibration 
data provided either on CCT-ADDS tapes In the Scrounge-era or on CCT-CALDUMP 
tapes In the TIPS-era. These studies characterize the In-orbit sensor 
performance, prior to ground processing. Scrounge-era digital products are 
described (P:Barker and Gunther). 

Background on radiometric procedures used In the Scrounge-era 
processing of the reflective bands Is provided (P:Barker, Abrams; (procedure 
paper)). Processing of Landsat-4 data during Scrounge-era was done one scene 
at a time. Results for approximately fifty Landsat-4 scenes are suimnarized In 
terms of wlthln-scene types of systematic variability and between-scene 
apparent changes In gain or sensitivity (S: Barker, Vol. I, 140-180, and 
P: Barker). Emphasis In these studies Is placed on characterizing precision, 
or relative radlometry, rather than on absolute accuracy. The most 
significant types of wlthln-scene systematic varlablly on the order of 1-2 DN, 
are: 

e Bin-Radiance Dependence 

e Scan -Corrlea ted Shifts 

e Coherent Noise 

§ Bright-Target Saturation. 

There may also be some variability, less than 0.5 DN, from: 

e Wlthln-Llne Droop 

e Forward/Reverse Scan Differences. 

The following papers deal with the radiometric characterization of 
In-orbit TM sensor performance from raw data: 

e P: Barker and Abrams, Characterization Paper 

Stability of Internal Calibration System 
Pre-Launch Sensor Performance (Vacuum Shift) 

- Post-Launch Sensor Performance 
Stability with Time 
Noise (Coherent) 

Location of Calibration Collect Window 
Non-Uniformity In the A/D Converter (Bin-Radiance 
Dependence) 

Recommendations for Operating and Processing 



• P: Lansl ng 

Thermal Band Characterization 

• P:Murphy 

Raw Data from the Canada Centre fo Remote Sensing (CCRS) 
Gains and Offsets 

Histograms (Bin-Radiance Dependence) 

Striping 

• I: Murphy 

Background In Shutter Region (Scan-Correlated Shift) 

Internal Calibration Region 

• PrJustlce 

Band-to-Band Correlations 

§ P:Metzler and I:Ma111a 

Scan Angle Effects (Wlthln-Llne Droop) 

Quantization-Level Histograms (Bin-Radiance Dependence) 
Down-Track Trace of kan-Llne Mean Signal (Scan-Correlated 
Shift) 

Correlation of Scan-Line Means 

• P:K1effer 

Variation of Response with Scan Direction (Forward/Reverse 
Scan Differences) 

Level Shift (Scan-Correlated Shift) 

Overshoot 
Gain Sag 

High Frequency Noise (Coherent Noise) 

t P: Barker ("R and R Paper") 

Table of Contents, as listed In Table 2* 

The characterization of possible wlthln-llne droop and 
foward-reverse-scan differences Is currently limited by an Inadequate 
quantification of the effects of bright target saturation (P:Barker). 

It Is assumed that between-scene variability of detector sensitivity 
or electronics will be calibrated out by adjustment In gain and offset to the 
raw data. One method to check the reasonableness of this assumption Is to 
monitor the apparent change In gain, and offset, with time and If It Is 
changing significantly, to see If that change can be accounted for by expected 
changes In a channel within a band, rather than a change In characteristics of 
the Internal calibration system. Results for Landsat-4 TM/PF apparent gain 
changes with time are presented and analyzed In the "R and R" paper (P:Barker). 

Image Data Quality 

If the sensor and calibration were perfect and unaffected by 
geometric resampling, then a user of TM digital Imagery need only acquire a 
calibrated and geometrically corrected CCT-PT tape and apply It directly to 
their Information extraction requirements. That Is the Ideal objective for 
the product of ground processing. This section Identifies the radiometric 
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Investigations of the final CCT-PT, and more generally, Its Immediate 
pre-cursor, the radlometrlcally unresampled CCT-AT. Both of these tapes are 
available to the public. To the extent that these Investigations of Image 
data quality find the radlometry within required tolerances, there Is no need 
to proceed backwards to determine whether or not the observed radiometric 
variability, or uncertainty, can be reduced by alternative ground processing 
procedures, or calibration, or whether It Is an uncorrectable. Innate 
characteristic of the sensor. Within the specifications of the system. It was 
the design objective of TIPS to produce a CCT-PT that would require no further 
calibration by the user. Specifications call for a post-calibration 
radiometric precision of ♦ 1 ON. In the absence of both a line-by-line 
calibration and an as yet undemonstrated capability for wl thin-line removal of 
effects of bright- target saturation, It will not be possible to obtain the 
objective In all scenes. Suggestions have been made for possible correction 
of known sources of systematic errors. 

The following papers report on the radiometric characteristics of 
radlometrlcally calibrated TM digital Image products: 

e P: Irons 

- TM Data Processing (During Scrounge-Era) 

e PrLyons 

Scrounge-Era Processl ng 

e I:F1sche1 

TM Radiometric Correction Algorithms (During Scrounge-Era) 

e P: Barker and Abrams, Procedure Paper 

Ground Processing (During Scrounge-Era) 

e P:Barker and Gunther 

- TM Image Products (Format and Characteristics of Tapes) 

e P: Murphy and I:Murphy 

Radiometric Processing at CCRS 

e P: Bernstein and I: Bernstein 

Sensor Data Entropy 
HI stograms 

A/D Converter Non-Linearity (Bin-Radiance Dependence) 

Failed Detectors (Landsat-) and Compensation 
Sensor Noise (32 KHz) and Reduction 
Radiometric Correction Processing 
Striping Removal 
Probabilistic Calibration 

• P:Bartolucc1 

Classification 

Clustering 

- Principal Components 
Covariance 

Correlation Matrices 
Water Temperature Mapping 



• I : Anuta 

- Coherent Noise (3.12 and 17 pixel/cycle) 

- Scan Angle Response 

• Principal Components 

- Temperature Mappl ng 

e P:Metz1er 

Qauntlzatlon Level Histograms (Bin-Radiance Dependence) 
Down-Track Means (Scan-Correlated Shifts) 

e P:K1effer 

- Detector- to-Detector Variations (Striping) 

- Scan-to-Scan Variations (Banding) 

Forward/Reverse Scan Differences 
Level Shifts (Scan-Correlated Shifts) 

Dynamic Response 

Gain Sag 

- Coherent Noise 
Between-Band Correlation 

e PrPodwysockl 

- Striping# 


SECTION II: GEOMETRY 

The TM geometric analyses/dlscusslons can be divided Into four 
categories: (1) pre-launch characterization of the sensor geometry. (2) 
discussions of the geometrical processing algorithms, (3) post-launch 
characterization of the geometrical behavior of the TM Instrument and (4) 
post-launch geometrical characterization of fully corrected (P-type) TM data. 
Under the first category, the geometry and geometrical measures of the TM 
sensor are discussed. The TM geometric correction operations and the 
differences between the Interim TM processing system (Scrounge) and the final 
system (TIPS) are presented In the second section. Under the third category, 
the emphasis Is on characterizing the magnitude of satellite vibrations 
(jitter) and the modulation transfer functions. In the fourth category, 
analyses have concentrated on evaluating band-to-band registration accuracy 
and the Internal geometric consistency of the TM data. The referencing 
convention followed by Barker In the first section of this Introduction will 
also be followed here. 

Stringent specifications on the radiometric and geometric performance 
of the TM Instrument necessitated greater complexity In the Instrument than In 
previous MSS Instruments. Key new TM geometric features Include (S:Engel, 

Vol. I, 41-61): (1) two focal planes, one containing the silicon detectors of 

bands 1-4 and the second containing the cooled detectors (InSb and HgCdTe) of 
bands 5-7, (2) nominally 30 meter IFOV's In bands 1-5, 7 with 16 detectors per 
bands, (3) bi-directional cross-tracking scanning with a scan line corrector 
to maintain a perpendicular scan, and (4) an angular displacement sensor to 
record vibrations at the sensor. A summary of the pre-launch geometrical data 
collected on the TM Instrument Is presented (S:Enge1, Vol I, 65-89): (1) 
measured Instantaneous fields of view (line spread function half widths), 
without electronics effects, (2) focal plane step response characteristics 



(rise-time, overshoot, settling times), (3) squire wove responses, (4) 
band-to-band registration, and (5) mirror scan profiles. These pre-launch 
measurements have been used (Schueler, 1983: Markham, 1984) to provide 
estimates of the line spread functions and transfer functions for the Thematic 
Mapper. Along-scan estimates of the line spread function half-widths were 
approximately 36 meters for bands 1-5, 7 and 137 meters for band 6. 

The TM aeometrlc ground processing was Initially performed by the 
Scrounge system (P:Lyon) and transferred to the operational Thematic Mapper 
Image Processing System (TIPS) (P: Beyer) In the Fall of 1983. Both of these 
systems use spacecraft attitude, ephemerls, and mirror scan correction data 
transmitted from the satellite, combined with pre-launch sensor data (e.g., 
band-to-band spacing and detector electronic delays) to generate systematic 
correction data (SCO). In addition In TIPS, control point processing Is being 
used to remove bias and drift errors In the SCO, generating Geodetic 
Correction Data (GCD). Then, using the GCD or SCO the data are geometrically 
corrected, registered to a map projection and resampled by a cubic 
Interpolation algorithm to 28.5 by 23.5 meter pixels. The standard cubic 
convolution weights were used In the Scrounge resampling process until i April 
1983. New weights, as recommended by Park, 1983 which result In a smaller 
overshoot near edges In Images were then Implemented ( I:F1 schel ) . The TIPS 
system will also use the revised weights. 

Due to the pre-launch concern with high frequency vibrations of the 
spacecraft, especially due to the scan mirrors, the post-launch TM Instrument 
geometrical analyses were concentrated on jitter. Most of the observed jitter 
was at 7KHz In the a1ong««can direction and can be easily corrected. The 
along-track effect of this jitter Is to produce underlap or overlap between 
adjacent scans. Analyses conducted (S:Kogut, Vol II., 54; P:Kogut) Indicated 
that the range of under! ap and overlap was In the neighborhood of 0.7 pixels. 

Several Investigators are examining the spatial resolution 
characteristics of the TM P-data (P:Schowengerdt, P:McG111em). Although this 
does not correspond to the spatial resolution of the TM Instrument due to 
resampling and atmospheric degradations. It provides an overall system 
measure. These analyses are concentrating on estimating the sensor Line 
Spread Functions (LSF), or modulation transfer functions (MTF) based on the 
rendition of step-like or llne-llko features on the Earth's surface that 
Intersect the scan pattern at oblique angles. Early results pointed to a 
system resolution (LSF half-power points) of 39 meters In a Webster Co., Iowa 
scene (P:McG111em) which Is consistent with sensor characteristics and a 
moderate degradation due to atmospheric and resampling effects. Late** 
estimates Indicated LSF half-widths of 39-45 meters (I: Anuta). The forward 
and reverse scan MTF's have also been compared as well as the MTF of the 
Interpolation/ resampling processing (I:Wr1gley). Results Indicated little. 

If any, difference In the system MTF from forward to reverse and a resampling 
MTF that was as expected. 

Band-to-band registration accuracy was the most Intensively studied 

? some trie characteristic of TM P-data (P:Card, PrFusco, P:Gurney, P:Podwysdkl, 
:Yao, IrBemsteln, I : Walker). Registration between the four primary focaf 
plane bands (1-4) was reported to be well within specifications (< .2 pixels) 
as was the registration between bands 5 and 7. Between the prime focal plane 
and the secondary focal plane, a misregistration of 0.5 to 0.8 pixels along 
scan and 0.2 to 0.3 pixels across-scan was reported. In addition, band 6 data 
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was found to be mis registered by approximately 3 resampled pixels. Shifts In 
the locations of the focal planes due to launch stress and vibration testing 
(bands 1-4 versus 5 and 7) and software errors (band 6) were deemed 
responsible for the misregistrations. Software changes were Incorporated to 
account for the shifts and errors In early 1983 (I:F1schel). 

The lack of geodetlcally corrected data has limited the ability of 
Investigators to characterize the geodetic accuracy of processed TM P-data. 
Early results found that the discrepancy between systematically corrected data 
and the ground Is largely translational (S:Thormodsgard, Yol. II, 62-63) 
Indicating that the Internal geometry Is quite good and that the Images can be 
corrected with only a low-order mapping to a projection. Later results 
^:We1ch) have confirmed the high Internal accuracy of TM data. 

Additional geometrically related analyses of TM data Included (1) ai 
analysis of the geometric accuracy of the GSFC fllmwrlter used for TM Images, 
which showed that TM low frequency distortions are likely to be less than the 
distortions Introduced by the fllmwrlter (P:Batson) and (2) a linear feature 
detection analysis which confirmed the resampled pixel size as 28.5 meters ard 
showed linear feature detection capabilities to as small as 6 meters 
(P:Gumey). 

Overall, the TM Instrument and ground processing techniques appear to 
be functioning properly and to within specification to the extent tested. 
Pre-launch concerns of excessive jitter did not materialize. Post-launch 
estimates of the system's spatial resolution are consistent with pre-launch 
estimates. In the Initial Scrounge tapes there were band-to-band 
ml sregl strati on si^he source of these errors has been located and apparently 
corrected. The scrounge system did not geodetlcally correct the data, 
however, systematically corrected data Is Internally consistent and the 
discrepancy with ground projections appears to be primarily translational 
and/or rotational. The TIPS processing system currently In use Is under 
evaluation for Its geometric properties: band-to-band registration and 
geodetic accuracy In particular. 
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